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Abstract 

Oh! In this article we consider Levy driven continuous time moving average pro- 

cesses observed on a lattice, which are stationary time series. We show asymptotic 
normality of the sample mean, the sample autocovariances and the sample autocor- 
relations. A comparison with the classical setting of discrete moving average time 
series shows that in the last case a correction term should be added to the classical 
Bartlett formula that yields the asymptotic variance. An application to the asymp- 

• totic normality of the estimator of the Hurst exponent of fractional Levy processes 

O ' 

m 



is also deduced from these results. 
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1 Introduction 



Statistical models are often written in a continuous time setting for theoretical 
reasons (e.g. diffusions). But if one wants to estimate the parameters of these models, 
one usually assumes only the observation of a discrete sample. At this point a very 
general question, the answer of which depends on the model chosen, is to know if the 
estimation should not have been performed with an underlying discrete model in 
the beginning. In this article we will consider this for moving average processes and 
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we refer to the classical moving average time series models as a discrete counterpart 
of this continuous model. 

To be more specific, let L = (L t )tem be a two sided one-dimensional Levy process, 
i.e. a stochastic process with independent and stationary increments, cadlag sample 
paths and which satisfies Lq = 0. Assume further that L has finite variance and 
expectation zero, and let / : R — > R be in L 2 (R). Let fj, € R. Then the process 
(X t ) t £R, given by 

X t = fi+ [ f(t-s)dL s , t€R, (1.1) 

can be defined in the L 2 sense and is called a continuous time moving average process 
with mean fi and kernel function f, driven by L. See also [6] for more information 
on such processes, in particular fractional Levy processes. The process (X t ) te R is 
then strictly stationary. Equation (jl.ip is the natural continuous time analogue of 
discrete time moving average processes 

X t = fx + J2^t-iZi, t€Z, (1.2) 

where (Z t )t<=z is an independent and identically distributed (i.i.d.) noise sequence 
with finite variance and expectation zero, and (ipi)iez is a square summable sequence 
of real coefficients. The asymptotic behaviour of the sample mean and sample au- 
tocorrelation function of Xt in (jl.2p has been studied for various cases of noise 
sequences (Zi)i<=z, such as regularly varying noise (cf. Davis and Mikosch [7]), mar- 
tingale difference sequences (cf. Hannan [9]), or i.i.d. sequences with finite fourth 
moment or finite variance but more restrictive conditions on the decay of the se- 
quence (cf. Section 7 of Brockwell and Davis [3]). 

Another approach to obtain limit theorems for sample autocovariances is to 
prove strong mixing properties of the time series under consideration, and provided 
it has finite (4 + <5)-moment, use the corresponding central limit theorems (such 
as in Ibragimov and Linnik |10j . Theorem 18.5.3). If even stronger strong mixing 
conditions hold, then existence of a fourth moment may be enough. Observe however 
that processes with long memory are often not strongly mixing, and in this paper 
we are aiming also at applications with respect to the fractional Levy noise, which 
is not strongly mixing. 

In this paper we shall study the asymptotic behaviour as n — > oo of the sample 
mean 

n 

Xn-A —n-^XiA, (1.3) 
of the process (X t )t£M. defined in (jl.ip when sampled at (An) ng N, where A > is 



2 



fixed, and of its sample autocovariance and sample autocorrelation function 

n—h 

% A (Ah) := n" 1 ^(A iA -X nA )(X (i+h)A -X nA ), h G {0, . . . , n - 1}, (1.4) 

i=l 

p n - A (Ah) : = 7n;A(A/l)/7 n;A (0), h G {0, . . . , n — 1}. (1.5) 

We write N = {0, 1, 2, . . .}. Under appropriate conditions on / and L, in particular 
assuming L to have finite fourth moment for the sample autocorrelation functions, it 
will be shown that X n . A and (/5 n;A (A), . . . , p nA (hA)) are asymptotically normal for 
each h G N as n — > oo. This is similar to the case of discrete time moving average 
processes of the form f)1.2|) with i.i.d. noise, but unlike for those, the asymptotic 
variance of the sample autocorrelations of model (jl.ip will turn out to be given by 
Bartlett's formula plus an extra term which depends explicitly on the fourth moment 
of L, and in general this extra term does not vanish. This also shows that the "naive" 
approach of trying to write the sampled process (X nA ) n£ z as a discrete time moving 
average process as in (|1.2p with i.i.d. noise does not work in general, since for such 
processes the asymptotic variance would be given by Bartlett's formula only. If 
p = 0, then further natural estimators of the autocovariance and autocorrelation 
are given by 

n 

j* n;A (Ah) := n-^XiAX^A, fc€{0,...,n-l}, (1.6) 

i=l 

p* nA (Ah) := 7* nA (Ah)/ 7 * n . A (0), fc€{0,...,n-l}, (1.7) 

and the conditions we have to impose to get asymptotic normality of 7*. A and p* n . A 
are less restrictive than those for "f nA and p nA . 

We will be particularly interested in the case when / decays like a polynomial, 
which is e.g. the case for fractional Levy noises. For a given Levy process with 
expectation zero and finite variance, and a parameter d G (0,1/2), the (moving 
average) fractional Levy process (Mf.^teR with Hurst parameter H := d + 1/2 is 
given by 



Mi 



t:d 



1 



(t-s) d + -(-s) d + 



dL„ tel (li 



r(d + i) 

(cf. Marquardt [II]). A process also called fractional Levy process was introduced 
before by Benassi et al. [2], where (x) + = max(x, 0) is replaced by an absolute value 
in dESJ, 

poo r 

dL s , t G R. (1.9) 



ML 



\t-s\ d - \s\ d 



Although both processes have different distributions, they enjoy similar properties. 
For instance the sample paths of both versions are Holder continuous, have the same 
pointwise Holder exponent, and they are both locally self-similar (see [2] for the 
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definition of this local property of their distributions). The corresponding fractional 
Levy noises based on increments of length A > are given by 

X\ = M\. d - Mt_ A . d , teM i = 1, 2. 

Hence the fractional Levy noise is a Levy driven moving average process with kernel 
function 

/l,A( s ) = f(^lj (s d + - (s - A) d + ) , s£l, (1.10) 

or 

fl A (s) = \s\ d -\s-A\ d , sem. (l.n) 

While the kernel functions f dA , i G {L2}, do not satisfy the assumptions we will 
impose for the theorems regarding the sample mean X n;A and the sample auto- 
correlation function p n - A , for d G (0,1/4) they do satisfy the assumptions we im- 
pose for the asymptotic behaviour of so that an asymptotically normal es- 
timator of the autocorrelation and hence of the Hurst index can be obtained if 
al £ (0, 1/4). For general d G (0, 1/2), one may take the differenced fractional Levy 
noises M\. d — 2M t ! _ A . rf + M^_ 2A . d , t G K, and our theorems give asymptotically 
normal estimators for the autocorrelation function of these processes. Please note 
that asymptotically normal estimators of the Hurst exponent for M?. d are already 
described in [2] but they use fill-in observations of the sample paths X 2 (k/2 n ) for 
k = 1, . . . , 2 n — 1. If L is a Brownian motion, then X 1 = CX 2 , where = means 
equality in distribution for processes and C is a constant, is the fractional Brown- 
ian motion and it is self-similar. Except in this case, fractional Levy processes are 
not self-similar and therefore observations on a grid k/2 n do not yield the same 
information as the time series X l (t), t G Z. 

The paper is organised as follows: in the next section we will derive asymptotic 
normality of the sample mean. Then, in Section [3] we will derive central limit the- 
orems for the sample autocovariance 7 n; A and the sample autocorrelation p n; A) as 
well as for the related estimators 7*. A and p* n . A of (jl.6p and (jl.7p . As a byproduct of 
the asymptotic normality, these quantities are consistent estimators of the autoco- 
variance and autocorrelation. Section [4] presents an application of our results to the 
estimation of the parameters of fractional Levy noises, where the underlying Hurst 
parameter is estimated. We also recall there that fractional Levy noises are mixing 
in the ergodic-theoretic sense, and we prove that they fail to be strongly mixing. 

Throughout the paper, unless indicated otherwise, L will be a Levy process with 
mean zero and finite variance a 2 = EL 2 , and X = (Xt)tzs. denotes the process 
defined in (jl.ip with kernel / G L 2 (M), / : K — > K. Its autocovariance at lag /i£K 
will be denoted by 

1 (h)= V (h) = Cov(X ,X h ) = a 2 [ f(- s )f(h-s)ds, (1.12) 

JR 
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where the last equation follows from the ltd isometry. 
Let us set some notations used in the sequel. 

If v is vector or A a matrix the transposed is denoted by v', respectively by A'. 
Convergence in distribution is denoted by -4. 

The function 1a for a set A is one for x E A, and vanishing elsewhere. 

The autocorrelation of X at lag h will be denoted by p(h) = pf(h) = 7(/&)/7(0)- 



2 Asymptotic normality of the sample mean 

The sample mean X n -/\ of the moving average process X of (jl.ip behaves like the 
sample mean of a discrete time moving average process with i.i.d. noise, in the sense 
that it is asymptotically normal with variance a 2 Yll^L-oo 7(^A), provided the latter 
is absolutely summable. 

Theorem 2.1. Let L have zero mean and variance a 2 , let p E R and A > 0. 
Suppose that 

oo \ 

F A : [0,A] -> [0,oo], u^F A (u)= £ |/(u + jA)| Ul 2 ([0,A]). (2.1) 

j=-oo y 

T/ien E^-oo l7(Aj)| < oo, 

2 



oo „a / 

/=-oo u \j 



£ /(u + jA) d«, (2.2) 



j=— oo \J =— °o 

and the sample mean of X A , . . . ,X nA is asymptotically normal as n —> oo, more 
precisely 

yfnX n - A A N I p, a 2 I I /(« + Aj) I du I as n — )• oo. 



Remark 2.2. Throughout the paper the assumption that L has zero mean can be 
dropped very often. For instance, if f G L 1 (M) n L 2 (R), i/ien the assumption of zero 
mean of L presents no restriction, for in that case L' t := L t — tE(L\), t G R, defines 
another Levy process with mean zero and the same variance, and it holds 

X t = p + E{L l ) [ f(s)ds+ I f(t-s)dL' s , t€Z, (2.3) 
which has mean p + E(L\) L /(s) ds. 

Proof. For simplicity in notation, assume that A = 1, and write F = F\. Continue 
F periodically on R by setting 

oo 

F(u)= \f( u + 3)l U ^ R - 

j=—oo 
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Since 

■2 



h f (h)\<a 2 \f(-s)\\f(h-s)\da 

J — oo 

by (fLT2jh we have 

-| OO n QQ OO 



h=—oo 00 h=—oo 

00 

|/( S )|F( a )da 

-00 

00 „^ 



£ / \f(s + j)\F(s)ds 
[ -F(s) F(s) ds < 00. (2.4) 

JO 



The same calculation without the modulus gives (|2.2p . 

The proof for asymptotic normality is now much in the same spirit as for discrete 
time moving average processes, by reducing the problem to m-dependent sequences 
first and then applying an appropriate variant of Slutsky's theorem. By subtracting 
the mean we may assume without loss of generality that \i = 0. For m G N, let 
fm ■= f l(-m,m), and denote 

Xi m) := / f m {s) dL s = [ t+m f(t -s)dL s , t G Z. 

Observe that (xj )t e z is a (2m— Independent sequence, i.e. (Xj )j<t and (^cj m ' > )j>t+2r 
are independent for each t € Z. From the central limit theorem for strictly station- 
ary (2m — Independent sequences (cf. Theorem 6.4.2 in Brockwell and Davis [3]) 
we then obtain that 

n 

=n- 1 / 2 J2 X t m) ^ Y(m) , n^oo, ( 2 -5) 
t=\ 

where y( m ) is a random variable such that 

Y^^N(0,v m ) 

with v m = E|=-2m7/ m (i)- Since lim m ^oo7/ m 0') = 7/0') for ea ch j € Z by (|1.12|), 
since 

,2 



l7/i ( _ m , m) (i)l<^ / I/O" 

and X^o \f(~ s ) \ I/O — s )l < 00 by (|2.4p . it follows from Lebesgue's domi- 

nated convergence theorem that linim^oo v m = X^-oo7/0)- Hence by (|2.2p . 



y (m) 4 Y, m ^ 00, where Y = N 0, cr 



jT 1 ( f; /(n + j)j duj . (2.6) 
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A similar argument gives lim^oo J^'jL-oo 7/-/ m (j) = °> so tnat 
lim lim Var fn 1/2 (X n -i -J&j) 

m— >oo n— >oo \ ' / 

( n />oo 
n / {f{t — s) — f m (t — s))dL s 

t=l J -°° 

oo 

lim V 7/-/ m (i) = 0, 



m— >-oo 

j=-oa 



where we used Theorem 7.1.1 in Brockwell and Davis [3] for the second equality. 
An application of Chebychef's inequality then shows that 

lim lim sup P(n 1/2 \X n;1 -X%i | > e) = 



for every e > 0. Together with (|2.5p and (12.6h this implies the claim by a variant of 
Slutsky's theorem (cf. [3], Proposition 6.3.9). □ 

Remark 2.3. Let us start with an easy remark on a necessary condition on the 
kernel f to apply the previous theorem. Obviously Fa G L 1 ([0, A]) is equivalent to 
f G L 1 (R). Hence F A G L 2 ([0, A]) =>■ / G L 1 (M). 

Remark 2.4. Unlike for the discrete time moving average process of (jl.2p . where 
absolute summability of the autocovariance function is guaranteed by absolute summa- 
bility of the coefficient sequence, for the continuous time series model (II. ID i£ is no£ 
enough to assume that the kernel satisfies f G L 1 (IR) n F 2 (M). An example is given 
by taking A = 1 and 

10, u < 0, 

1, UG[0,1), 
l:i " 2 V}' '' (» ./)'• ue[j,j + l), j g N. 
For i/ien i/ie function F\ is given by 

F 1 (u) = ^2f(u + j) = (l-u)-^ 2 , ue [0,1), 

so t/iflt Fi G ^([0,1]) \ L 2 ([0, 1]). 5u< Fi G ^([0,1]) is equivalent to f G L 1 (R), 
and since < 1 /or a// this implies also f G F 2 (IR). Observe further that 

for non-negative f , condition (|2.1|) is indeed necessary and sufficient for absolute 
summability of the autocovariance function. 
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3 Asymptotic normality of the sample auto- 
covariance 

As usual, we consider the stationary process 

/oo 
f(t-s)dL s , tel. (3.1) 
-oo 

We recall that 

n 
t=l 

and first we establish an asymptotic result for Cov (7*. A (pA), 7*. A (gA)). 

Proposition 3.1. Let L be a (non-zero) Levy process, with expectation zero, and 
finite fourth moment, and denote a 2 := EL\ and n := a' ~ A EL\. Let A > 0, and 
suppose further that f <E L 2 (R) n L 4 (R) and that 



l[0,A]->R, Y, f( u + kA f G^ 2 ([0,A]). (3.2) 

\ fc=— oo / 

For q G Z denote 

oo 

0, ;A : [0, A] -> R, ^ /(u + kA)f(u + (k + q)A), 

k=— oo 

which belongs to L 2 ([0, A]), fry ^ e previous assumption. If further 

oo 

£ | 7 (/*A)| 2 <oc, (3.3) 

fe=— oo 

f/ien we /iou e /or each p, q G N 

/•A 

lim nCov (7*. A (pA), 7« ; a(<?A)) = (»/ - 3)^ 4 / g P ;A(u)g q . A (u) du+ 

oo 

2 [7(fcA) 7 ((A;-p + g )A) + 7((fc + g )A)7((fc-p)A)]. (3.4) 

fc=— oo 

Proof. For simplicity in notation we assume that A = 1. The general case can be 
proved analogously or reduced to the case A = 1 by a simple time change. We shall 
first show that for t,p,h,q G Z 

/oo 
/(«)/(« + p)f(u + h+ p)f(u + h+p + q)du 
-oo 

+7(P)7(?) + 7(^ + P) 7 (fr + ?) + 7(^ + P + ?) 7 (>0- (3-5) 
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To show this, assume first that / is of the form 

m/e 

f{s)=fmA s )= ^* 1 (ie,(i+l)e](*)> ( 3 - 6 ) 

i=— m/e 

where m G N, e > such that 1/e 6 N, and ^ 6 I, i = —m/e, . . . , m/e. Denote 

poo m / e 
Xf, m , e ■■= J fm,e(t - s) dL s = ^ ^i(L t -ie ~ L t _( i+V)e ), t G R. 

%=— m/e 

Denote further 

Zi := Li e — L(i_i) e , ieZ. 
Then (Zi)i £ z is i.i.d. and we have 

m/e 

^te;m,e = 4>i Z t-i, t £ Z. 

i=—m/e 

At this point we will need to compute the fourth moment of integrals of the Levy 
process. Let us state an elementary result that yields a formula for this moment. 

Lemma 3.2. Let <fi £ L 2 (M) n L 4 (M) ; then, with the assumptions and notations on 
L used in Propositio ns. 11 



E([ (j){s)dL s f = (rj - 3)a 4 f cj) 4 (s)ds + 3a 4 ( f ct> 2 {s)ds) 2 . (3.7) 
Jr Jr Jr 

Proof. If v is the Levy measure of L and A its Gaussian variance, then by the Levy 

Khintchine formula we get 

£(it) = E exp ^iu J (f>(s)dL s ^j 

= ex V (-]-Au 2 f tf(s)ds+ [ [e iu ^ s)x - I - iu(f){s)x\v{dx)ds\ . 
V 2 7r JrxR J 

Then E( L cp(s)dL s ) 4: is obtained as the fourth derivative of £ at it = 0. If we recall 
that (77 — 3)<7 4 = f R x 4: u(dx), and a 2 = .A + f R x 2 is(dx), we get (13. 7j) . after elementary 
but tedious computations. □ 

To continue with the proof of Proposition 13.11 we now apply (|3.7p to the special 
case where f(s) = l( , € ](s) and we get 

EZ 2 = EL 2 = a 2 e, EZf = EL\ = r/a^e — 3cr 4 e + 3cr 4 e 2 . (3.8) 

As shown in the proof of Proposition 7.3.1 in [3], we then have 

^(-^t;m,e-^t+p;rn,e-^t+h+p;m,e-^i+h+p+g;m,e) 

m/e 

= [EZf - 3(EZ 2 ) 2 ) AA+p/elpi+h/e+p/elPi+h/e+p/e+q/e 
i=—m/e 

+lm,e{p)lm,e{q) + J m ,e( h + P)lm,e{h + q) + -f m ,e( h + p + q)^/ m>e (h), 
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where 7 m , e (u) = E(XQ. m ^X u . m ^), u G R. By ([HTH]), 

£Z 4 - 3(£Z 2 ) 2 = fa - 3)a 4 e, 

and 

m/e 

e l J J i' l Pi+p/e' l Pi+h/€+p/e' l Pi+h/e+p/iE+q/e 
i=—m/e 

X) 

f(u)f(u + p)f(u + h + p)f(u + h + p + q) du, 

oo 

so that (Eg) follows for / of the form / = f m>e . Now let / G L 2 (R) n L 4 (M) and X u 
t £ R, defined by (|3.ip . Then there is a sequence of functions (fm k ,e k )keN of the form 
(|3.6p such that fm k ,e k converges to / both in L 2 (R) and in L 4 (R) as k — > oo. Then for 
each fixed t £ R, we have that Xt- mk ^ k — > Xt in L 2 (P) (P the underlying probability 
measure) as k — > oo, where we used the ltd isometry. Further, by Lemma 13.21 and 
convergence of fm k ,e k both in L 2 (R) and in L 4 (R), we get convergence of X t;Tnkj e k to 
X t in L 4 (P). This then shows (|3.5p . by letting fm k ,e k converge to / both in L 2 (R) 
and L 4 (R) and observing that ^m k ,e k { u ) l(u) for each u £ R. From ()3.5[) we 
conclude that, with p,g£N, 

Cov( 7 ; ;1 (p), 7 ; i (g)) = n- 1 ^2(l-n- l \k\)T k , (3.9) 

|fc[<n 



where 



T fc = ^(k)^{k-p + q)+^(k + q)^(k-p) 

/•oo 

+ (r, - 3)<r 4 / /(«)/(« + p)/(« + k)f{u + q + k) du. 



Now by O, Er=~oo < 00 if 



E 



fc=— oo 

is finite. Denote 



f(u)f(u + + k)f(u + q + k)du 



(3.10) 



G P («) := \f(u + k)f(u + k + r)\, u£R, r G N. 

A;=— oo 

Then G r is periodic, and by assumption, G r restricted to [0, 1] is square integrable. 
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Hence we can estimate (|3.10p by 

x /'OO 

E / |/(«)/(« + + *:)/(« + ? + A;)|d« 

= E / |/H/(« + p)|Gf,(«)d« 

h=— 00 
00 „i 

= E / |/(« + /»)/(« + P + /i)|G ff (t*)d« 

= / Gp{u)G q {u) du < 00. 

The same calculation without the modulus and an application of the dominated 
convergence theorem to ([3.90 then shows f|3.4j) . □ 

Remark 3.3. A sufficient condition for (|3.3p is that Y2hL-oo < °°? which 

is implied by the function in Theorem \2.1\ belonging to L 2 ([0, A]). Another suf- 
ficient condition is that $ : u 1— > J] fcgZ |J r (/)(u + 27r/c/A)| 2 , n £ [0, 27r/A] is in 
L°°([0, 2-7r/A]), where J-(f) is the Fourier transform of / £ L 2 (R) in i/ie /orm 
z H> Jf^ e lzt f(t) dt (for L l -functions). For if ||$||oo < -B, f/ien (/(• + hA)) he z is a 
Bessel sequence in L 2 (R) with bound B/A, i.e. 

E / <p(u)f(u + hA)du < BA- 1 / (£(n) 2 dn V 93 £ L 2 (R), 
j, J —00 J— 00 

ra=— 00 

see e.g. Theorem 7.2.3 in Christensen J3]/. Taking <p = f then gives the square 
summability of the autocovariance functions by (|1.12p . 

Please remark that X^/j^=-oo 7(^^) 2 < 00 cannot be deduced from the condition 
that u 1 y ^Ca^=-oo/( u + ^A) 2 *s *n ^ 2 ([0i^])- One can take A = 1 and f(s) = 
Ei>i / or 5 < -ff < |, to i/ie /atter condition but not 7(^) 2 < 00 • 

Remark 3.4. 5y (|1.12p . i/ie condition (|3,3p can ^ e written as 
X^fcL-oo ( f( s )f( s + &A) ds^ < 00. TTie assumption (|3.1ip nserf in Theorem \3.5\ 
below is slightly stronger than (13. 3D . 6ui equivalent to (13. 3D i/ / > 0. 

The following theorem gives asymptotic normality of the sample autocovariance 
and sample autocorrelation and the related estimators 7*. A and /0*. A . 

Theorem 3.5. (a) Suppose the assumptions of Proposition \3A\ are satisfied and 
suppose further that 

2 



E U \mf(s + kA)\ds\ <oo. (3.11) 



k=— oo 

T/ien we Ziawe /or eac/i /i £ N 

v^(7; ; a(0) - 7(0), • • • ,1* n - A (hA) - 7 (/i))' 4 N(0, V), n -> 00, (3.12) 
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where V = (v pq ) P7q=0 ,...,h G R^ 1 '^ 1 is the covariance matrix defined by 



f A 

v pq = (V ~ 3 )°" 4 / 9p-A( u )SqA( u ) du+ 
Jo 

oo 

J2 [7(feA)7((fc-p + g)A)+7((fc + g)A) 7 ((fc-p)A)]. (3.13) 

k=— oo 

In addition to the assumptions of (a), assume that the function 

oo 

\f( u + 3 A )\ 

j=~oo 

is in L 2 ([0, A]). Denote by 
n-j 

7n;A(jA) = n' 1 ^2(X tA - X n;A )(X {t+j)A - X n . A ), j = 0, 1, . . . , n - 1, 
t=\ 



the sample autocovariance, as defined in (|1.4p . Then we have for each h £ N 

V^(7« ; a(0) - 7(0), • • • ,%A(hA) - j(h))' A N(Q, V), oo, 
where V = (y pq )p t q=a,...,h is defined by (|3,13p . 

(c) For j € N let p* n . A (jA) = 7^aO'A)/7^-a(0) a ^ = 7n ; A(jA)/7 n;A (0), 

i/te latter being the sample autocorrelation at lag jA. Suppose that f is not almost 
everywhere equal to zero. Then, under the assumptions of (a), we have for each 
h £~N, that 

V^(p* A (A)-p(A),...,p* A (hA)-p(hA))' AN(0,W), n^oo, (3.14) 



_ /-A 

Wij- A = w ijlA ^ / (gi. A (u) - p(iA)g . A (u))(g j;A (u) - p{jA)g . A (u)) du, 



where W = W A = (miy ; A)i,j=i,...,ft w given 

7-3) 
7(0) £ ,o 
and 

oo 

t%;A = E {p{(k + i)A)p((k+ j)A)+ p{(k-i)A)p{(k+ j)A) + 2p{iA)p{jA)p(kAf 

k=—oo 

-2p(iA)p(kA)p((k + j) A) - 2p(jA)p(kA)p((k + i)A)) 



£ (p((fc + OA) + - £)A) - 2p(zA)p(A;A)) 



it=i 



+ j) A) + - j) A) - 2p(iA)p(fcA)) 

is given by Bartlett's formula. If additionally the function u h-» 2~^j*L-oo l/( n + J A.) | 
is in L 2 ([0, A]), then it also holds that 

V^(pn-A(A)-p(A),...,p n . A (hA)-p(hA))' 4JV(0,W), n ^ oo. (3.15) 
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Proof. For simplicity in notation we assume again A = 1 in this proof. 

(a) Using Proposition 13.11 it follows as in the proof of Proposition 7.3.2 in [3], 
that the claim is true if / has additionally compact support. For general / and 
m G N let f m := fl 

(— m,m)- Hence we have that 
™ 1/2 (7: ;(m) (0) - 7m(0), • • . ,l* n . (m) (h) - lm {h))' A Y m , m -> oo, 

where j m is the autocovariance function of the process X t - m = f m {t — s) dL s , 

7n-r m )(p) = n_1 H™=i ^t;m^*+|);m the corresponding autocovariance estimate, and 

Y m = N(0,V m ) with V m = {v pq;m ) Piq=0 ,...,h and 

[-1 00 
v pq ;m = (v-3)o- 4 gp;( m )(u)g q -( m )(u)du+ ^ [lm(k)j m (k-p+q)+'y m (k+q)j m (k-p 

Here, 5 p; ( m )(«) = Efct-oo + k)f m (u + k +p), u G [0, 1]. 
Next, we want to show that lim m _ i , 00 = y. Observe first that 

oo oo 

9 P ;{m){u) = ^2 fm(u+k)f m (u+k+p) ->• ^ f (u+k) f (u+k+p) = g Ptl (u) =: g p (u) 

k=— oo fc=— oo 

almost surely in the variable m as m -> oo by Lebesgue's dominated convergence 
theorem, since -u i— )■ X^fcl-oo + + ^ + is m 

L 2 ([0, 1]) by (1321) and hence 

is almost surely finite. Further we have 

00 

\g P] (rn){u)\ < \f(u + k)f( U + k+p)\ 

k=—oo 

uniformly in u and m, so that again by the dominated convergence theorem we have 
that flp;( m ) — > g P in L 2 ([0, 1]) as m — > 00. Next, observe that 

/oo 
\f(s)f(s + k)\ds VmeN VfceZ. 
-00 

Since lim m _>. 00 j m (k) = j(k) for every k G Z, it follows from the dominated con- 
vergence theorem and fj3. 1 1 f) that (jm(k))kez converges in Z 2 (Z) to (7(&))fceZ- This 
together with the convergence of <7 p; ( m ) gives the desired linim-^ = V, so that 

Y m -4 Y, m — >■ 00, 
where Y = N(0, V). Finally, that 

lim limsupP(n 1/2 |7*. (m) (p)- 7m (p)-7*(p)+7(p)| > e) = V e > 0, p G {0, . . . , h} 

follows as in Equation (7.3.9) in [3]. An application of a variant of Slutsky's theorem 
(cf. [5], Proposition 6.3.9) then gives the claim. 

(b) This follows as in the proof of Proposition 7.3.4 in [3]. One only has to observe 
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that by Theorem l2.ll ^JnX n -\ converges in distribution to a normal random variable 
as n — > oo. In particular, X Tl .\ must converge to in probability as n — > oo. 
(c) The limit theorem follows as in the proof of Theorem 7.2.1 in [3], and for Wij we 
have the representation 

Wij- a = (vij - p{i)v 0j - p{j)vio + p(«)p(j>oo)/7(0) 2 
(V ~ 3)^ 4 

= W a-A+ 7(0)2 x 
i 

{9i( u )9j{u) - p(i)g (u)gj{u) - p(j)gi(u)g (u) + p{i)p{j)go{uf) du, 
giving the claim. □ 



o 



Remark 3.6. It is easy to check that Wij-A = S^j;A if f is of the form f = 
S£-oo ^il(iA,(i+l)A] j i- n accordance with Bartlett's formula, since then (X tl \) t <=z 
has a discrete time moving average representation with i.i.d. coefficients. 

Remark 3.7. Another case when Wij- t A = Wij-& is when r\ = 3, which happens if and 
only if L is Brownian motion. However, in general we do not have Wij-^ = Wij-&. 
An example is given by f = l(o,i/2] + 1(1,2] an d A = 1, in which case g\-i = l(o 1/2] 
and go-i = 2 • 1(0,1/2] + 1(1/2,1] > an d it is easy to see that gi-\ — p(l)go-i is not almost 
everywhere zero, so that uin-i 7^ w\i-i if tj 7^ 3. The latter example corresponds to 
a moving average process, which is varying at the scale \, but sampled at integer 
times. Observe however that ^n ; i/2 = ^ii ; i/2 by Remark \3.b\ A more detailed study 
of such phenomena in discrete time can be found in Niebuhr and Kreiss [13]. 



Remark 3.8. Recently, sophisticated and powerful results on the normal approx- 
imation of Poisson functionals using Malliavan calculus have been obtained. E.g., 
Peccati and Taqqu 1151 Theorems 2, 3 and 5] prove a central limit theorem for 
double Poisson integrals and apply this to a specific quadratic functional of a Levy 
driven Ornstein-Uhlenbeck process, and Peccati et al. \14\ Section 4] obtain bounds 
for such limit theorems, to name just of few of some recent publications on this 
subject. It may be possible to apply the results of \14[ to obtain another proof 
of Theorem 13.51 under certain conditions such as finite 6th moment, but we have 
not investigated this issue further. Note that our proof uses only basic knowledge of 
stochastic integrals and methods from time series analysis. 



4 An application to fractional Levy noise 



We will now apply the previous results to fractional Levy processes. Recall from 
(jl.8p and (jl.9p that these were denoted by 

1 



r(d + i) 



(t - s)% - {-s)% 



dL s , t G 



and 
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ML 



oo 

d \ n \d 



dL=, t e 



|t — s\ — \s\ 

' — oo L 

respectively, and the corresponding fractional Levy noises based on increments of 
length A > by 

Xi = M}. d -Ml A . d , f £l, i = l, 2. 

Hence the fractional Levy noises are Levy driven moving average processes with 
kernel functions 



flA(s) = f ^ rr) (si-(s-A)i), seR, 

and 

fl A (s) = \s\ d -\s-A\ d , sem, 

respectively. Neither f dA nor /J A are in L 1 (M), so Theorem 12.11 cannot be applied 
because of Remark 12.31 Please note that for the same reason the assumptions for 
(b) of Theorem 13.51 and for (|3. 15j) are not fulfilled. 

For simplicity in notation we assume A = 1, and drop the subindex A. Although 
the fractional noises X % have different distributions for i = 1 and i = 2, they are 
both stationary with the autocovariance 

nxUhXi) = i*{h) = (\h+ - +\ h - ii m+i ) , (4.i) 

where C l (d) is a normalising multiplicative constant depending on al. Both processes 
X 1 are infinitely divisible and of moving average type, hence we know from [H [8] 
that (XDt^z is mixing in the ergodic-theoretic sense. For fixed h G Z, define the 
function 

F : R z -)■ M, {x n ) n& ' ^ xox/,. 
If T denotes the forward shift operator, then 

F(T (Xl) teZ ) = X l k X l k+h , 

and from Birkhoff's ergodic theorem (e.g. Ash and Gardner [TJ, Theorems 3.3.6 and 
3.3.10) we know that 

^Y j Xixi +h ^E(F{{Xi) t& ))=EXiXl n^oo, 
k=i 

for i = 1, 2, and the convergence is almost sure and in L 1 ([TJ, Theorems 3.3.6 and 
3.3.7). Hence, with 7* = 7*^ as defined in fjl .6[) . 

lim 7 * (/l) = C |fc + i|2<i+i _ 2 \h\ 2d+1 + \h- l\ 2d+1 ) a.s. 
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Since 7(0) = C*(cf)<7 2 and 7(1) = C i {d)a 2 {2 2d - 1), p* n {l) = is a strongly 



is a strongly consistent estimator for d. 

The question of the asymptotic normality of these estimators arises naturally. 
There are many classical techniques to show the asymptotic normality of an ergodic 
stationary sequence by assuming some stronger mixing assumption. As far as we 
know, they do not work in our setting. To illustrate this point, we would like to 
show that fractional Levy noises are not strongly mixing. Let us first recall the 
definition. 

Definition 4.1. Let (X n ) ne z be a stationary sequence, and let 

a x (n) = sup{|P(vl D B) - P(A)P(B)\, A G a(X k , k < 0), B G a(X k , k > n)}. 
The sequence (X n ) n< zi is strongly mixing if lim n _ i . 00 ax(n) = 0. 

In our case we know the weak mixing property linin^oo \P(A(~]B) — P(A)P(B)\ = 
for A G a(Xo), B G a(X n ), because of [H HJ- There are classical central limit 
theorems for strongly mixing sequences, see e.g. [T2] for an overview. The following 
result, which is stated as Proposition 34 in |12j . will be particularly useful for us. 

Theorem 4.2. Suppose that (Xt)t£Z is a mean zero, strongly mixing sequence and 
that there exists some 5 > and a constant K > such that 



consistent estimator for 2 — 1. Hence 





E\Xq 



2+6 



< 



00 



(4.3) 




V m G N. 



(4.4) 



(4.5) 




S^(t) :=Y,^ neN, fG[0,l] 



where [x\ is the integer part of the real number x. Then 




weakly in D[0, 1] 



(4.6) 



where B is a standard Brownian motion. 
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Lemma 4.3. Let L be a two sided non-zero Levy process with expectation zero and 
finite fourth moment, and let d G (0, 1/2). Let 

Xj:=Mi. d -MU;d, t€Z, 

be the corresponding fractional Levy noises, for i = 1, 2. Then (Xl) te % satisfies 
(Ojl - (03]) with 5 = 2, fori = 1, 2. 



Proof. The proof is only written for the fractional noise X 1 denoted by X but it is 
similar for X 2 . Equation (|4.3p holds since L has finite fourth moment and since the 
kernel function of fractional noise is in L 2 (R) n L 4 (M), and fl33| follows from the 
fact that 

m 

Y J X i = M m , d , (4.7) 
i=i 

and 

Var (M m . d ) = Cm 2d+1 VmeN (4.8) 
for some constant C. To see (j4.5|) for (5 = 2, we use (|4.7p and 

fm{s) := T ^ +l) [{m-s) d + -{-s)% S GK. 

Then by Lemma 13.21 

E\M d {m)\ 4 = (r]-3)a A [ fi{s)ds + 3a A ( [ fi(s)ds 



Observe that 



and that 



2 

r> , a j \ _ Cm 4d+2 



(Jj 2 m (s)ds^' 
J 



f^s)ds<C'm id+1 

for positive constants C, C, which gives the claim. □ 

A consequence of the previous theorem and lemma is the following negative 
result. 

Corollary 4.4. Assume the assumptions of Lemma \4-3\ Then the fractional Levy 
noises X 1 and X 2 are not strongly mixing. 



Proof. If fractional Levy noises were strongly mixing, then (|4.6p would follow. Please 
remark that, since fractional noises are increments, S^ n \t) = M d {\nt\) — Af^(0), 
and (Var (^™ =1 Xj)) = Cn 2d+1 by (|4.8p . Owing to the asymptotic self-similarity of 
fractional Levy processes (Proposition 3.1 in [2]), we know that 

M 2 ([n t\) - M 2 (0) d 

* B d + l/2{t), 



n" 
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where the limit is a fractional Brownian motion with Hurst exponent d + 1/2. A 
similar asymptotic self-similarity holds for M}. Hence (|4.6p is violated and the 
corollary is proved by contradiction. □ 

Nevertheless one can apply Theorem 13.51 to get the asymptotic normality of the 
estimator (14. 2|) . 

Proposition 4.5. Assume the assumptions of Lemma \4-S\ and let d be defined by 
(|4.2p . If d £ (0,1 /4) then \/n{d — d) converges in distribution to a Gaussian random 
variable as n — >■ oo. 

Proof. The proof is only written for the fractional noise X 1 denoted by X but it is 
similar for X 2 . 

We shall apply (|3.14p to get convergence of \/n{p^ l (V) — p(0)) to a Gaussian 
random variable. First, observe that 

fl, eL 2 (R)nL A (M), (4.9) 

and YlkL-oo 7(^) 2 < 00 • The latter inequality is classical for the fractional Gaussian 
noise, when d < 1/4, and holds for fractional Levy noises since they have the same 
autocorrelation as fractional Gaussian noise. This also implies (|3.1ip by Remark l3.4l 
since > 0. Let us check that go := go-i G L 2 (0, 1). Since 

oo 

T 2 (d + l)g (u) = Yl ((« + k )% ~ ( u + k ~ X )+) 2 ' 

k=— oo 

it follows for all u G (0, 1) that 

oo 

r 2 (d + i) 50 ( u ) < £((i + *)}. - (* - 

k=0 




< OO, 



so that even go £ L°°([0, 1]). Hence the assumptions of Theorem [33] (a) are fulfilled, 
and the result follows from (|3.14p . □ 

If d > 1/4 then ]CfcL-oo7(^) 2 = oo, since fractional Levy noises have the same 
autocorrelation as the fractional Gaussian noise. Hence we consider 

zi = xi-XU, tez, 

for which it holds ^2 k L_ 00 Jz i (^) 2 < 00 • We conclude from Birkhoff's ergodic theo- 
rem that 

1 n 

7nz(h) ■= ~ y^Z k Z k+h -> E(Z Z h ), n-^oo, 

k=2 
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is a strongly consistent estimator of 

E(Z Z h ) = (\ h + 2 \ 2d+1 + A\h + l\ 2d+l - 6\h\ 2d+1 +±\h- l| 2d+1 -\h- 2\ 2d+ ~ 



Therefore p* n z (l) = ^"' z | Q j is a strongly consistent estimator for 4>{d) = 3 + s ^ 2 +2 + — -• 
It turns out that 4> is increasing on (0, 1/2). Therefore we can define the estimator 

d:=r\p*n,z^))- (4-10) 

Proposition 4.6. Assume the assumptions of Lemma \4-S\ and let d be defined by 
(|4.1Up . If d G (0, 1/2) f/ien yjn(d—d) converges in distribution to a Gaussian random 
variable as n — >■ oo. 

Proof. The proof is only written for the fractional noise X 1 denoted by X but it is 
similar for X 2 . Let us remark that Z t = ~~ s ) ^si where 

/i>) = /i>)-/i>-i), 

To apply Theorem 13.51 we have to check that 

fU s )eL 2 (R)nL*< 



which is obvious from (|4,9p . Moreover we already know that Ylh=-oo 7z i {k) 2 < oo. 
This time, however, the kernel function is not nonnegative, but it is easy to see 
that 1/^0)1 < Cmin(l, \t\ d ~ 2 ) and hence that 

/oo 
\fll(f) fd,l(t + k)\dt<C min(l, {hi* 1 - 1 ), VfcGZ, 
-oo 

for some constants C,C, giving (|3. 1 1[) . Finally, 

oo 

T 2 (d + l)g (u) = (( n + k )i + ( u + k ~ 2 )+ _ 2 (n + k — l) d + ) 2 , 

k=—oo 

and estimating the summands separately for k < 0, k = 0, 1 and k > 2 we obtain 
for u € [0, 1] 



T 2 (d+l)g (u) < l + (2 d + 2) 2 + Y,{k d -(k-l 

< oo, 



d\2 

n- - V a, - i; J 

fc=2 



so that <7o £ ^°°([0, 1]) C L 2 ([0, 1]). The claim now follows from Theorem 13. 5[ using 
(I3T41) . □ 
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